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ABSTRACT 

Hydrogen fluoride has been established to be an excellent tracer of molecular 
hydrogen in diffuse clouds. In denser environments, however, the HF abundance 
has been shown to be approximately two orders of magnitude lower. We present 
Herschel/HIFI observations of HF J = 1 — toward two high-mass star formation 
sites, NGC6334 I and AFGL 2591. In NGC6334 I the HF line is seen in absorption 
in foreground clouds and the source itself, while in AFGL 2591 HF is partially in 
emission. We find an HF abundance with respect to H 2 of 1.5 • 10" 8 in the diffuse 
foreground clouds, whereas in the denser parts of NGC6334 I, we derive a lower 
limit on the HF abundance of 5 • 1CT 10 . Lower HF abundances in dense clouds 
are most likely caused by freeze out of HF molecules onto dust grains in high- 
density gas. In AFGL 2591, the view of the hot core is obstructed by absorption 
in the massive outflow, in which HF is also very abundant (3.6 ■ 1CT 8 ) due to the 
desorption by sputtering. These observations provide further evidence that the 
chemistry of interstellar fluorine is controlled by freeze out onto gas grains. 

Subject headings: stars:formation — ISM: molecules 
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Introduction 



Hydrogen fluoride (HF) is an exceptional molecule, because of its peculiar chemistry. 
Fluorine is one of the few atoms which has a greater affinity to hydrogen than hydrogen 
itself, and thus the reaction 

H 2 + F -> HF + H 

is highly exothermic (AE w 16000 K, although the activation of this reaction energy is 
400 K and quantum mechanical tunneling through this barrier has to be considered). The 
destruction of HF occurs by the relatively slow photo-dissociation (1.17 • 10~ 10 s _1 ) and 
reactions with ions, such as He + , Hg , and C + , which are rare compared to H 2 . Thus HF 
formation is much more efficient than its destruction under most interstellar conditions. As 
a consequence HF is considered to be the main reservoir of fluorine in molecular regions, 
i.e, where H2/H > 1, and therefore it is thought to be an ideal tracer of H2. Astrochemical 



models predict an abundance of HF relative to H2 (a 



paper are relative to H 2 ) of 3.6 -10 8 in such clouds (INeufeld et al 



1 fractional abunda nces given in this 



20051 ). 



Neufeld et al. 



(119971 ). who observed the 



The first detection of HF was reported by 
HF J = 2 - 1 transition (v = 2463.43 GHz) with relatively low resolution (R=9600) 
in absorption towards Sgr B2, using the Infrared Space Observatory (ISO). However, 
because high densities or strong radiation fields are needed to populate the HF J = 1 
level, the HF column density derived from this observation is subj ect to large uncertainties. 



With the advent of the HIFI instrumen t (jde Graauw et al. 



Space Observatory (IPilbratt et al. 



20101 ) aboard the Herschel 



20101 ) . high resolution spectroscopy (R > 10 6 ) of 
the ground state rotational transition of HF has become possible for the first time. 
Several HF observations in diffuse cloud along the lines of sight towa rd strong sub- 



millimeter dust continuum sources, such as W31C (INeufeld et al. 



(jSonnentrucker et al 



20101), and Sgr B2(M) (IMonie et al. 



2010), W49N and W51 



201 ll ) have been reported. 
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HF J = 1 — is detected in absorption toward all of these continuum sources, which was 
expected considering the extreme conditions necessary to populate the excited rotational 
states of this molecule. These observations prove that HF is an excellent tracer of H 2 
in diffuse clouds, although the abundances of HF relative to H2, derived by comparison 
with CH, are ~ 1.5 ■ 1CT 8 , with an uncertainty of 46%, and thus a factor of two below 
the chemical model predictions. Furthermore, HF turned out to be a much more sensitive 
tracer of H2 than the commonly used CO, which is in agreement with the model predictio n 



that H 



2005|). 



is more abundant t 



Sonnentrucker et al. 



ran CO in diffuse clouds with low extinction (INeufeld et al 



(j2010l ). for example, detected a diffuse cloud on the line of sight 



toward W51, not seen in CO. 



The only HF abundance calculated in a high -mass star-forming r egion, and not in 



Phillips et al. 



fl2010h . who detected 



a diffuse foreground cloud, has been reported by 
HF J = 1 — in absorption towards Orion KL. Orion KL is a peculiar source, heated from 
the front, which shows hardly any lines in absorption at submillimeter wavelengths. The HF 
abundance derived by Phillips et al. is a lower limit of 1.6 • 10~ 10 , approximately two orders 
of magnitude lower than the HF abundances reported in diffuse clouds. This result can be 
either explained by geometrical effects, due to the incomplete coverage of the continuum 
source by absorbing material, or, under the assumption that most fluorine is bound in HF, 
is caused by HF molecules sticking onto the surfaces of dust grains under high-density 
conditions. The only three det ections of HF in emis sion so far are in the immediate vicin ity 



of th e AGB-star IRC+10216 (jAgundez et al. 



20111). the Orion Bar ([Van der Tak et al. 



20121 ) and towards the Seyfert 1 Galaxy Mrk 231 (IVan der Werf et al 



2010) 



To investigate the abundances of HF in dense regions in more detail, we analyze here 
the HF J = 1 — spectra towards two high-mass sta r-forming reg ions, NGC6334 I and 



AFGL 2591. NGC6334 is a relatively nearby (1.7 kpc; 



Neckel 



19781 ) high-mass star-forming 



-6 - 



region, which harbors sites of many stages of protostellar evolution (IStraw & Hyland 



19891 ). Single-dish continuum observations at su b-millimeter w avelength show that a total 



mass of 200 M & is associated with NGC6334 I fjSandell 



NGC6334 I, studie d extensively over the last decades (e.g. 



2000 ). The molecular hot core 



Beuther et al 



Hunter et al. 



20061 ). is associated wit h an ultra compact HII region (jde P 



and shows a very line-rich spectrum (jSchilke et al. 



2006 



2008 



2007; 



ree et 



Thorwirth et al. 



observations in NGC6334 I have led to the detection of CH (IVan der Wie 



(jEmprechtinger et al. 



20101 ). H 2 0+ flOssenkopf et al.l 120101 ). and H 2 C1+ flLis et al 



et al 



al. 



19951 ) 



2003 ). HI FI 



20liH 2 O 



201 



showing that 



lydrides are common tow ard this source. Furthermore, a bipolar outflow 



(Leurini et al. 2006; 


Beuther et al. 


2008 


) and 


H 9 0, OH, CH,OH class II, and NH, masers 


have been detectec 




Kraemer & Jackson 


1995: Ellinesen et al. 1996: Walsh et al. 


2007). Rolffs et al. 


(2011 


studied the radial structure of several high-mass star-forming 



cores, including NGC6334 I. This investigation, mainly based on HCN observations with 
the APEX telescope, revealed a density law of n oc r~ L5 . The temperature also follows 
a power law, with an index of 0.83 and 0.4 in the inner and the outer part, respectively. 
Interferometric data (SMA), however, revealed that the internal structure of NGC6334 I 
is much more complex. The hot core consists of four compact condensations within a 10" 
region, emitting ab out 50% of the con tinuum flux, whereas the other 50% stem from the 



extended envelope (jHunter et al. 



2006|). 



AFGL 2591 is an isolated high-mass star-forming region located within the Cygnus-X 



region. The distance toward AFGL 2591 was reported to be betwee n 0.5 kpc and 2.0 



and for modeling purposes a distance of 1 kpc was assumed (e.g. 



kpc, 



Van der Tak et al. 



19991 ). However, new measurement base on VLBI parall ax measurement s of H 2 masers 



20121 ). Based on a 



suggest that the distance might be as large as 3.3 kpc (IRygl et al. 
distance of 1 kpc, the luminosity of AFGL 2591 is ~ 2 • 10 4 L Q , however, adopting the 
newly measured distances the luminosity would be approximately a factor ten higher. The 
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embedded central star, which is obscured in the visible by the mas sive envelope, has an 



estim ated mass of 16 M and an effective temperature of 33,000 K (IVan der Tak &: Menten 



20051 ). Besides the circumstellar envelope a massive outflow with an outflow velocity o f 



20 kms 1 with respect to the systemi c velocity has been detected (ILada et al. 



In addition, 



19841 ). 



Hasegawa fc Mitchelll (1l995[ ) found a second, faster (v out ~ 40km s 1 ), but 



weaker outflow. Many attempts to model the structure of AFGL 2591 have been carried 
out a ssuming spherical symm e try and a power-law density and tempera t ure structure 



Van der Tak et al 



1999 



Van der Tak et al 



2000 



Dotv et al 



200 9]) . A cavity, caused by the outflow might by present fjVan der Tak et al. 



2002 



de Wit et al. 



Van der Wiel et al. 



( 120111 ) mapped AFGL 2591 in the frequency range from 330-373 GHz 
with the JCMT. In this frequency range they found ~ 160 lines, of which 35 show extended 
emission (> 15"). These observations reveal a small scale structure (< 10 4 AU), and a line 
of sight velocity gradient is apparent in most molecules. Their modeling suggests that a 
non isotropic structure or a velocity gradient must be present on ~ 10 4 AU scales in order 



to explain the observations. 



Observations 



The observations, used in th is paper, are part of the Chemical iffirschel Surveys 



of Star forming regions (CHESS; 



Ceccarelli et al. 



20101 ). a Herschel guaranteed time key 



program, which has conducted unbiased spectral line surveys of several star-forming 
regions. The HF J = 1 — spectra presented in this paper are observed in two 
sources: NGC6334 I (« 200 o : 17 /l 20 m 53.32 s , /3 200 o : -35°46'58.5") and AFGL 2591 
(a 2 ooo : 20 h 29 m 24.9 s , /3 20 oo : +40°11'21.6"). The observation identification (ObsID) 
number of the scans are 1342206594 and 1342196510 for the observations in NGC6334 I 
and AFGL 2591, respectively. Data were taken on May 12 th 2010 (AFGL 2591) and 
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October U th 2010 (NGC6334 I). 



Both observations were carried out with the HIFI instrumen t (Ide Graauw et al. 



on board of the Herschel Space Observatory ( iPilbratt et al. 



2010|) 



20101 ) using the double b eam 



switc h mode (180" chopper throw), and have been reduced with the HIPE pipeline ( jOtt 



20101 ) version 5.1 and version 3.0, respectively. The double sideband (DSB) spectra were 
observed with a redundancy of eight, which allows the deconvolution and isolation of the 



single sideband (SSB) spectra ( jComito fc Schilke 



20021 ). In the case of NGC6334 I the 



deconvolution was done in HIPE and the deconvolved SSB spectra were exported to the 
FITS format for subsequent analyses in the IRAM/GILDAS package Q. The spectra of 
AFGL 2591 where exported to the FITS format right after the basic data reduction (level 
2), and further data reduction and deconvolution was done using GILDAS. The spectra 
shown in this paper are equally weighted averages of the H and V polarization. 

The HF J = 1—0 line {y = 1232.47627 GHz) has been clearly detected in both sources. 
The beam size of Herschel at this frequency is 18" and a main beam efficiency of 0.62 was 
applied. Both spectra are displayed in Figure [TJ 



3. Results & Analysis 

As shown in Figure [IJ the HF line seen is in absorption, in both sources, over a 
wide range of velocities. Only in the velocity interval from —8 kms -1 to —2 kms -1 in 
AFGL 2591 a weak emission is seen. In many cases the HF absorption is completely 
saturated. Assuming that all HF molecules are in the rotational ground state, one can 



www.iram.fr/IRAMFR/GILDAS / 
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calculate from the absorption depth the line opacity and subsequently the HF column 
density using following formulae: 

Tub 



T 



-In 



-^~back , 

where T MB is the measured brightness temperature and T back the temperature of the 
background continuum, and 

37r 3 / 2 ■ Av gi 



(1) 



N. 



HF 



■ T 



(2) 



2Vln2 A 3 ■ A9u 

where A is the wavelength (243.2 /im), A the Einstein A coefficient (2.422 • 10~ 2 s" 1 ), and 
gi and g u are the statistical weights of the lower and upper state (1 and 3, respectively). 
The assumption of all HF molecules being in the ground state is justified by the hig h 



201l|) 



Einstein A coefficient and the low rates for collisional excitation (IGuillon et al. 
corresponding to a critical density of n cr i t = 1.9 • 10 8 cm" 3 , and the relatively high energy 
of the J = 1 state (50 K). As shown by Van der Tak (submitted), most alternative, 
non-thermal excitation mechanisms, such as infrared pumping and formation pumping, 
play only a negligible role. Excitation caused by submillimeter radiation on the other hand 
is possible in the hot inner part of high mass star forming regions and has to be considered. 
We will discuss the potential influence of this mechanism on the HF excitation individually 
for each source later on. 




v LSR (kms ) 



v lsr ( kms ') 



Fig. 1.— HF J = 1 - spectra in NGC6334 I (left) and AFGL 2591 (right). The respective 
continuum level as well as the zero level are marked with the gray lines. 
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3.1. NGC6334 I 



Many of the features seen in the HF spectrum toward NGC6334 I originate from four 
foreground clouds, of which two reside in the immediate environment of NGC6334 I, whereas 
the two remaining clouds are not associated wit h the source itself. Thes e foreground 
clouds, which have been identified in CH spectra (IVan der Wiel et al.ll2010l ) are at velocities 
of -3.0 kms -1 , 0.0 kms -1 , +6.5 kms -1 , and +8.0 kms -1 , whereas the velocities of the 
individual components of NGC6334 I itself are between —6.0 kms -1 and —8.0 kms -1 . 
From fitting several C 18 lines observed with HIFI, we know that the source itself can be 
modeled by three components at t>LSR=~6.0 kms -1 , -8.0 kms -1 , and -6.5 kms -1 (Plume 
priv. communication). The first two co mponents are ident ified as two of the embedded 



sub cores (SMA 1 and SMA 2 following 



Hunter et al. 



20061 ). whereas the later component 



represents the lower density envelope. Because, HF is a sensitive tracer of the molecular 
gas, we expect to see all these components in the HF spectru m as well. A massiv e outflow, 



detected in water (Emprechtinger et al. submitted) and CO (jLeurini et al. 



2006|) in the 



velocity range between -100 km sand +60 kms, is not seen in HF, neither in emission nor in 
absorption. The upper limit of the optical depth of HF in the blue lobe is 0.88 what results 
in an upper limit on the HF column density of 8.42 • 10 13 cm -2 and therefor e in an upper 



Leurini et al. 



l imit o n the HF abundance of 5 • 10 , using the H 2 column density given by 
(120061 ). The upper limit of the integrated intensity of the red lobe is 2.69 Kkms -1 , which 
results, assuming the canonical HF abundance, in an upper limit on the HF excitation 
temperature of 10 K. Both limits are consistent with standard HF abundances and negligible 
HF excitation. Thus the non-detection of the outflow in HF is another indication that most 
HF molecules are in the ground state in this source. The foreground components and the 
envelope are expected to be more extended than the continuum sources in NGC6334 I, 
and thus T back in equation [TJ is equal to the measured continuum level in the spectra. The 
components associated with embedded cores, however, absorb only continuum emitted 
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from the respective core. Furthermore, we cannot be certain that the absorbing material 
covers all the continuum radiation emitted by the core. Therefore, even at very high optical 
depths only a fraction of the background radiat ion is absorbed. From interferometric 



observations at A = 1.3 mm (IHunter et al. 



20061 ) we estimate the contribution of each of 



the cores to be approximately 25% of the total continuum flux. This 25% value is only 
a rough estimate, since the actual value depends on the exact dust temperature and the 
coverage of the continuum source. Therefore we use the actual continuum contribution 
of the cores as free fitting parameters f& and f C 2- Using the centroid velocities and line 
widths of the seven components derived from the CH and C 18 observations, and assuming 
that all HF molecules are in the ground state, the only free parameters in addition to 
f c i and f C 2 are the HF column densities of the individual components. This approach is 
valid for all components except the +8.0 kms -1 foreground cloud, which is narrower in HF 
(Av = 2.5 kms -1 ) than in CH (Av = 3.5 kms -1 ). Using equations [I] & [2] with the free 
parameter f c i, f& and the HF column densities as model to fit the observed spectra, we 
used the X 2 method to determine the column densities listed in Table [TJ results are shown 
in Figure [2J The reduced X 2 for the best fit result is 1.07, but all results with a reduced X 2 
lower than four have to be considered as possible solutions. Thus the HF column densities 
can be determined only with an accuracy of 20%-50% for each individual component. 



The HF column density of the kms - and the +8.0 kms - foreground components are 
well determined (±20%), whereas the —3.0 kms -1 component is less constrained (±50%) 
due to the blending with the absorption features from NGC6334 I itself. For the +6.5 kms -1 
component only a lower limit for the column density can be given since this line is optically 
thick. However, the opacity of the +6.5 kms -1 component is moderate, because we do not 
see any effects of extremely high optical depth, such as a line width wider than the optically 
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Table 1: Fit results for the individual velocity components in NGC6334 I. 



Component 


V 


Av 


N 




[kms -1 ] 


[kms -1 ] 


[cm -2 ] 


Foreground 


-3.0 


5.0 


2.4+1.2 ■ 10 13 


Foreground 


0.0 


3.0 


1.5 + 0.3 ■ 10 13 


Foreground 


+6.5 


1.2 


> 1.4 ■ 10 13 


Foreground 


+8.0 


2.5 


4.0 + 0.8 • 10 12 


Envelope Cloud 


-6.5 


4.0 


2.7+1.3 • 10 13 


SMA 1 


-6.0 


6.0 


2.0 + 1.0 • 10 14 


SMA 2 


-8.0 


5.0 


2.0 + 1.0 • 10 14 


f cl = 0.32 ± 0.06 




f c2 = 0.21 ± 0.06 





t 1 1 1 1 1 1 1 1 1 r 




J I I I I I I I I I L 

-20 20 

Velocity (km/s) 



Fig. 2. — Comparison of the normalized spectrum of HF (black) and the model (red). The 
green lines show the individual components. 
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thin CH line, or a saturated line shape. Two narrow residuals seen at vlsr = ~1 kms -1 and 
-9 kms -1 , are not reproduced by our fit. Considering that both residuals are located in 
the wings of major absorption features, a possible explanation is a velocity or line width 
gradient within individual components. These two features may also be caused by emissioji 
lines from other specie s. However, there a re no plausible candidates in the JPL database 



or the CDMS catalog flMuller et all 



20051 ) CH is known to be an excellent tr acer of H? at 



densi ties below 10 4 cm 3 with an abundance of 3.5 • 10 8 with respect 



Van der Wiel et al- 



t o H? ( jSheffgr et al. 



fl2010h to derive the 



20081 ). Using the CH column densities reported by 
H 2 column density, we are able to determine the HF abundance relative to H 2 in the diffuse 
foreground clouds (see Table [2]) . 



The HF abundance in the clouds at -3.0 kms 1 , 0.0 kms 1 , and +6.5 km s 



-i 



agree fairly well with the abundances found in other diffuse clouds (iMonje et al. 



Sonnentrucker et al. 



2011 



20101 ) However, the +8.0 kms 1 component exhibits an HF abundance 



which is significantly lower. The +8.0 kms 1 component is also the component in which 
the HF line is much narrower than CH. Furthermore, the abundance of CH + is significantly 



2 http:/ /spec.jpl. nasa.gov/home. html 



Table 2: Abundances of HF in the foreground components in NGC6334 I 



vel. Comp 


iV(HF) 


iV(H 2 ) 


X(HF) 


(kms -1 ) 


(cm -2 ) 


(cm -2 ) 




-3.0 


2.40 • 10 13 


1.8 • 10 21 


1.3 • 10 -8 


0.0 


1.50 • 10 13 


6.0 • 10 20 


2.5 ■ 10 -8 


+6.5 


> 1.4 ■ 10 13 


8.6 • 10 20 


> 1.6 ■ 10 -8 


+8.0 


4.0 ■ 10 12 


6.0 ■ 10 20 


6.7- 10 -9 
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enhanced in this cloud compared to the other foreground components (Lis et al. in prep, 
i ndicating; that the chemical and physical properties of the gas may be very different 



( IGodard et al 



20091 ) 



Since observations of Orion KL suggest that the abundan ce of HF in dense g as in 



20101 ). we 



high-mass star-forming regions is lower than in diffuse clouds (jPhillips et al. 
attempted to derive the HF abundance in the components associated with NGC6334 I 
itself. The fits for the SMA 2 core and the envelope material are relatively good and HF 
column densities can be determined with an accuracy of 50%. The column density of the 
SMA 1 core, however, is poorly constrained since this component is blended with several 
other components. As can be seen in Figure [2j the HF absorption features originating 
from SMA 1 and SMA 2 are largely overlapping, which results in a large ambiguity of 
the individual continuum coverage factors (±0.2). However, their sum can be defined 
quite well as 0.53 ± 0.05. The fit results for SMA 1 presented in Table [1] are based on 
the assumption that the HF abundance relative to H 2 in SMA 1 and SMA 2, which have 
similar physical conditions, is the same. To estimate the H 2 column density, we are using 
the C 18 column densities derived by Plume (priv. communication) from an LTE fit to 
seven C 18 transitions. We are assuming standard values for the C 16 0/C 18 ratio (500) 
and the C 16 abundance (9.5 ■ 10~ 5 ). Furthermore we divide the resulting column density 
by two, because approximately half of the material emitting the CO lines resides behind the 
continuum sources, and therefore does not contribute to the HF absorption. The resulting 
HF abundances range from 3.7 • 10~ 10 in the envelope to 5 — 8 • 10 -10 in the dense cores. 
These values are approximately two ord ers of magnitude lo wer than the abundances in 



diffuse clouds a nd predicted by mod els (INeufeld et al. 



limits found by 



Phillips et al. 



20051 ) but comparable to the lower 



( 120101 ) in Orion KL. In these calculations we did not consider 



submillimeter radiative excitation, and therefore these results are lower limits. However, 



taking the observed continuum values (IHunter et al 



20061 ) and using LVG calculations, we 
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found that we underestimate the HF column density by less than a factor four for conditions 
present in the dense cores. Therefore, the low abundances of HF derived in NGC6334 I 
itself cannot be explained by HF excitation alone. 

The calculated HF abundances are based on the assumption that approximately half of 
the mass associated with NGC6334 I resides in front of the corresponding continuum source. 
However, the actual structure of the source is not known, and thus the HF abundance might 
be underestimated. But, it should be noticed that it is very unlikely that the unknown 
geometry of NGC6334 I is the sole reason for the low HF abundance derived here. Assuming 
that only geometrical effects account for the low HF abundance would imply that in all 
three components (SMA 1, SMA 2, and the envelope cloud) as well as in Orion KL, the 
bulk of the molecular gas (~ 98%) resides behind the main continuum source. 



3.2. AFGL 2591 

The spectrum of AFGL 2591 is much simpler than that of NGC6334 I; in addition to 
one foreground absorption component at ~ kms" 1 , only one broad component ranging 
from —6 kms -1 to —20 kms -1 can be seen. In the velocity range between the two 
absorption features a weak HF emission is detected. Both absorption features are most 
likely very optically thick, since all background radiation is completely absorbed in the 
corresponding velocity ranges. Unfortunately, the continuum of AFGL 2591 is much weaker 
than the one of NGC6334 I and thus the signal to noise ratio is much lower, resulting in a 
lower limit for r of 1.6. 

The profile of the HF J = 1 — line is very peculiar and only the profiles of the 



water ground-state transitions p-H 2 In — Ooo and p-H 2 2 i2 — loi 
although the emission component is much stronger in the water lines 



ook somewhat similar 



Choi et al 



20111). The 
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broad absorption component is clearly blue shifted with respect to the systemic velocity of 
-5.7 kms -1 of AFGL 2591, and thus most likely caused by the massive outflow . This outflow 



is also seen in CO and HCO + (ILada et al 



1984 : 



Hasegawa &: Mitchell 



1995 



, and centers 



at a velocity of -7 kms -1 , slightly blueshifted with respect to the envelope (IVan der Wiel 
201ll ). He also derived an H2 column density of the outflow of 1.5 ■ 10 22 cm -2 , based on 
observations of several CO line and radiative transfer calculations. This H 2 column density 
leads to an HF column density of 2 ■ 10 14 cm -2 assuming an HF abundance of 1.5 • 10 -8 . 
Using a line width of 10 kms -1 , which is the approximate velocity width of the outflow, 
this HF column density corresponds to a r of ~ 9 (using equation [2]), which is in agreement 
with the high optical depth seen in the observed spectrum. HF abundances of the order 
of ~ 3 • 10 -10 , as found in star-forming regions, are clearly too low to cause the observed 
deep absorption. This is a very interesting aspect, because the temperature (100-150 K) 
and density (10 5 — 10 6 cm -3 ) are comparable to the conditions in high-mass star-forming 
regions. It should be noted that a shock with a velocity of ~ 15 kms -1 , as derived form the 
spectrum, may ca use a partial dissociation of H? (up to 80%) due to energetic collisions with 



neutral particles (IFlower &: Pineau Pes Forets 



20101 ) . This may cause an underestimation 



of the total hydrogen column density in the outflow (2N(H 2 )+N(H)). However, collisions 
with enough energy to destroy H 2 can potentially destroy HF and CO, too. Furthermore, 
in the absence of H 2 , the HF production is hindered. Therefore, we are able to derive the 
HF abundance relative to H 2 , but not relative to the total hydrogen density. The error of 
the HF abundance relative to H 2 is determined by the uncertainty of the H 2 column density 
derived from CO. We estimate that these errors are definitely less than the degree of the 
dissociation of H 2 and therefore clearly lower than one order of magnitude. 

The foreground absorption line width (equivalent width is 3.2 kms -1 ) is clearly l arger 



than the corresponding feature in the CO lines (Av = 0.8 kms ; 



Van der Wiel 



more comparable to the width of the feature seen in CH (Av = 1.8 kms l . 



201 lh. and 



Bruderer et al. 
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2 1 Oh . Using the line width of CH as the Av caused by Doppler broadening for HF, and 
fitting the absorption feature using equations (|T]) and (j2j), leads to a minimum optical depth 
of eight, and thus a minimum HF column density of 4 • 10 13 cm -2 in the fo reground clou d. 



Comp aring this HF column density with the H2 column density found by 



Van der Wiel 



( 1201 ll ) based on CO observations, yields an HF abundance relative to H2 of > 1.3 ■ 10~ 8 , 
which is a typical value found in diffuse clouds and close to that found in the clouds toward 
NGC6334 I. However, the C O abundance in diffu s e clou ds might be significantly lower than 



the canonical value of 10 4 . 



Sonnentrucker et al. 



( 120071 ) found CO abundances with respect 



to H 2 ranging from 2.75 ■ 10~ 7 to 2.75 ■ 10~ 5 in material on the line of sight toward several 
stars. Hence the total H2 column density toward AFGL 2591 may be underestimated 
and, because we derive only lower limits on the HF column density, the uncertainty of the 
derived HF abundance is very large. 

To get a better understanding of the HF spectrum in AFG L 2591, especially of the out 



flow a nd the HF emission, we model this source using RATRAN ( jHogerheijde &: Van der Tak 



20001 ). RATRAN is a Monte Carlo radiative transfer code, which includes excitation by 
collision of the molecules with H2 and excitation by submillimeter radiation. It calculates 
a grid of spectra integrated along a pencil beam, which are subsequently convolved to 
the desired spatial resolution. The Einstein coefficien ts of the modeled transitions are 



taken from the LAMDA database (jSchoier et al 



rates are used ( 



Guil 



Van der Wiel et al. 



on et al. 



20051 ). whereas newly calculated collision 



201lh . The source model is based on the "static model" of 



(120 111 ), who modeled the envelope of AFGL 2591 assuming a power law 



density distribution (re = r L0 ) and used a self-consistent solution for the dust temperature 



( [Van der Tak et al. 



20001). To take the outflow 



expansion velocity to each shell in the model of 
with radius as 



component into account , we add an 



Van der Wiel et al. 



( 1201 ll ). which increases 



J exp 



i-roy 

r / 
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where Voo is the maximum velocity of the outflow and ro is the radius at which the outflow 
starts. /3, which was set to two, is a measure of how quickly the maximum velocity 
is reached; r was set to a value so that the H 2 column den sity outside th is radius is 
comparable to 1.5 • 10 22 cm -2 , the N# 2 found in the outflow (IVan der Wiel 



201 lh . At radii 



smaller than ro the expansion velocity is zero. The maximum velocity is set to 9 kms -1 
to match the width of the absorption feature. The modeled HF spectrum is not very 
sensitive to the actual velocity profile, and a proper solution can be found for profiles with 
moderately increasing v exp . The foreground component is fitted using the foreground option 
in RATRAN assuming an optical depth of eight and a line width of 1.8 kms -1 . 

Deriving the HF abundance is somewhat difficult, because the absorption features 
are optically thick. Hence we calculate models with three different abundances relative to 
H2. In Model 1 we assume an abundance of 1.5 • 10~ 8 , the HF abundance found in diffuse 
clouds, whereas in Model 2 an HF abundance of 5 ■ 10~ 10 as found in the dense part of 
AFGL 2591 is adopted. In Model 3 we assume a value of 5 ■ 10 -10 in the inner part of 
the source, where v exp = kms -1 and 1.5 • 10~ 8 in the outer part. A comparison of the 
three models with observations is shown in Figure |3j Models 1 and 3 fit the observations 
quite well, whereas an HF abundance as low as 5 • 10 -10 in the outflow can be ruled out. 
The model calculations show as well that the emission seen at wlsr ~ — 4 kms" 1 , stems 
from the red lobe on the far side of the continuum source, and not from the hot core itself. 
The emission at vlsr ~ — 4 kms -1 is overestimated by both, Model 1 and Model 3. This 
may be due to geometry effects, since the outflow is modeled by a shell, whereas it is in 
fact a collimated outflow. The absorption of all radiation in the blue lobe of the outflow 
indicates that most HF molecules are in the J = state, and thus excitation by collisions 
with H2 or any other effective excitation mechanism can be ruled out. As suggested by our 
model calculation, the emission from the red lobe is most likely caused by absorption and 
immediate re-emission of 1.2 THz photons from the continuum source of AFGL 2591. 
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Fig. 3. — Comparison of the observed HF J = 1 — transition (black) with radiative 
transfer models. Model 1 (red): HF abundance is 1.5- 10 -8 ; Model 2 (green): HF abundance 
is 5 • 10~ 10 ; Model 3 (blue): HF abundance is 5 • 10~ 10 in the quiescent gas and 3.6 • 10~ 8 in 
the outflow. 
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Because none of the features of the HF spectra stems from the hot quiescent part of 
AFGL 2591, we cannot distinguish whether the abundance in the inner part is similar to 
the diffuse cloud abundance or to the lower limits fou nd in the high-mass star-forming 



region NGC6334 I and Orion KL ([Phillips et al. 



20101 ) . The assumption of a uniform, low 



HF abundance in the dense parts of high-mass star-forming regions is consistent with our 
data. 



4. Summary & Discussion 

We determine the abundance of HF in very different physical environments, and obtain 
very different results. In agreement with previous studies, we find the HF abundance of 
~ 1.5 ■ 10~ 8 in diffuse clouds, implying that about half of the available fluorine is in the form 
of HF in the gas phase. An exception is the +8.0 kins" 1 component toward NGC6334 I, in 
which HF is less than half as abundant as in other diffuse clouds. The chemical composition 
of this foreground component, however, seems to be different from the other components 
(e.g., strongly enhanced CH + abundance; Lis et al. in prep.), which requires further studies. 

In the high-mass star-forming regions, at much higher densities compared with diffuse 
clouds, the HF abundance seems to drop by about two orders of magnitude. The correlation 
with the physical conditions, especially density and temperature, would be crucial to 
determ ine the nature of this HF depletion, but it still remains to be found 



Phillips et al 



(120101 ) argues that the low abundance of HF at higher densities is due to freeze out of HF 
onto dust grains. This argument is drawn from the fact that almost all fluorine should 
be bound in HF due to its high proton affinity, and that the desorption energy of HF is 
assumed to be quite large due to its polar nature. In dense (10 5 — 10 6 cm -3 ) and warm 
(100-150 K), but dynamically active regions such as the outflow in AFGL 2591, we find 
HF abundances close to the values found in diffuse clouds. In the picture of HF freeze 
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out at high densities this can be explained by HF desorption du e to sputtering by high 



velocity particles, similar to the desorption of water in outflows f lKristensen et al 



Emprechtinger et al 



2010h 



2010 ; 



In both sources HF is the only fluorine bearing species we find in the HIFI line survey, 
and we can confirm non-detections for CF, CF + , DF, and H 2 F + . In NGC6334 I the upper 
limit for CF + J = 5 — 4, a mol ecule which has bee n detected in the Orion Bar with an 



abundance of a few times 10 10 (iNeufeld et al 



20061 ). is 0.012 K. Assuming cloud properties, 



as determined by C 18 fits (Plume priv. communication), yields an upper limit of the CF + 
abundance of 6 ■ 10 -12 . This low upper limit, 50 times lower than the abundance found in 
the Orion Bar, can be explained by the different nature of the two objects. C + , a precursor 
of CF + , and thus CF + itself, are expected to be much more abundant in photo-dominated 
regions, such as the Orion Bar, than in massive cores. The upper limit of the abundance 
of these four fluorine bearing species together is ~ 7 ■ 10 -11 , and thus these molecules hold 
less than 0.4% of the total fluorine. The non-detection of other fluorine bearing species is 
another indication that freeze out onto dust grains causes the low HF abundances observed 
in dense, quiescent gas. 

This analysis demonstrates that the utility of HF as a tracer of H 2 , as predicted by 
chemical models, is only applicable to diffuse clouds. The fluorine chemistry in dense clouds 
is not yet fully understood and further theoretical and observational studies are required. 
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